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' In this talk I discuss the non-linear development of magnetic fields in the early 

universe. Since this is based on a classical field theory, which turns out to have a 
^ ! rather complex structure, I thought it could be of interest for this meeting, under the 

H ' heading of "applied field theory" . I very briefiy mention a number of particle physics 

mechanisms for generating magnetic fields in the early universe, but the emphasis is on 
the field theoretic aspects of the developments of these fields, mainly on the occurence 
of inverse cascades, i.e. generation of order from disorder. In this connection it is also 
discussed how the Silk effect (photon diffusion) is counteracted by the inverse cascade, 
which moves energy from smaller to larger scales. 

Many galaxies (including our own) are observed to have magnetic fields. One way 
to observe such fields is to study the polarization of light passing the galaxy. Due to the 
interaction with the field and the plasma there is a Faraday rotation of the polarization 
vector, proportional to the field and to the square of the wave length of the light. In 
this way fields are found to have the order of magnitude 10~^ — 10^^ Gauss on a scale 
of 100 kpc|^ If you have forgotten what a G(auss) is: the mean field on the sun is 
approximately one G. 

Usually the galactic magnetic field is explained by the dynamo effect: turbulence 
(e.g. differential rotation) in the galaxy enhances the magnetic field exponentially up to 
some saturation value, corresponding to equipartition between kinetic and magnetic en- 
ergy. The dynamics which governs these phenomena is called niagnetohydrodynamics, 



"A p(arse)c is an astronomical unit, which has the physical value 1 pew 3.26 light year. 
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abbreviated as MHD, which is essentially the Maxwell plus Navier-Stokes equations. 
The dynamo can produce an enhancement factor of several orders of magnitude. An 
important feature is that the dynamo needs a seed field. It appears reasonable to 
assume that this field is of primordial origin, i.e. it has existed already in the early 
universe. Astrophysicists often say as a joke that a primordial magnetic field is a field 
which has existed for so long that everybody has forgotten how it was created. How- 
ever, in particle physics we must be more serious since we have knowledge of the early 
universe, and hence we should explain the origin of these fields. 



PRIMORDIAL MAGNETIC FIELDS IN THE EARLY UNIVERSE 

In natural units magnetic fields have dimension (mass)^. At the electroweak scale, 
assuming the Higgs mass to be of order mw, there is essentially only one mass, mw, 
and we may therefore expect something like 

Bew < ml, ^ 10^4 G (1) 

on a scale ~ l/rriw. This is a huge field, far larger than anything one has ever seen or 
produced on this earth. How does this compare with the rather weak fields found in 
galaxies? 

In the standard cosmological model all distances are blown up by the scale factor 
R{t). It is useful for estimates that the scale factor is proportional to the inverse 
temperature. Thus, Ruow/Rew = ^Ew/^now ~ 10^^. Hence, an initial correlation 
length of order ~ l/rriw is of order 1 cm today, which has no astrophysical interest. 
We need fields on a scale of order 100 kpc ~ 3x10^^ cm. 

If we assume that B is essentially random, we can estimate the field at any distance 
from a simple random walk. We have the field at the initial correlation length, but we 
want it at ~ 10^^ times this length. Thus, in d dimensions we have 

< Bew >scale 10-/n^v.^ G/(1023)^/l (2) 

So for (i = 3 we get < -Bew >~ 10^^°G, whereas for li = 2 and d = 1 we have < -Bew > 
approximately equal 10 G and lO^^G, respectively, on the scale of lO^^/mi^/. 

In order to see if these fields are reasonable, we need to know the cosmological 
developments of < -B >. From MHD (with viscosity ignored) one has the result that 
the flux through a surface bounded by a curve following the fluid of charged particles 
is conserved. Since such a surface increases like R^, it follows that < B > decreases 
like 1/-R^ when the universe expands with the scale factor R [|. It therefore follows 
that if today we need e.g. a primordial field of order 10~^^G on a scale of 100 kpc, 
then on the corresponding scale 10"^^ /mw at the electroweak phase transition, we need 
< B >^ 10^^ G. Thus, from the random walk estimates above we see that only the case 
d = 1 comes near this value, although a factor 10^ is missing. Actually one could argue 
that the case (i = 1 is relevant, because in observing the magnetic field by Faraday 
rotation, a one dimensional average is made along the line of sight. However, this 
argument is not really convincing, since the dynamo effect is three dimensional, and 
hence the field relevant for this effect is the very small 3d average. 



^The metric is 



1 — fcr^ 

where k = +1, 0, —1 for a closed, flat or open universe, respectively. 



(3) 
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The conclusion is thus that if fields of the order of can be generated at the 
EW-scale, then there could still be missing a factor of order 10"^, where x is of order 
3. However, it should be emphasized that a field of order m^r is very large, and is not 
obtained in most mechanisms for creation of primordial fields. Hence, in most cases x 
is larger than 3. 

MAGNETIC FIELDS FROM PARTICLE PHYSICS MECHANISMS 

In this section we very briefiy discuss a number of proposed mechanisms for the 
generation of primordial fields. The list is by no means exhaustive. 

Fields From Inflation 

An inflationary creation of primordial magnetic flelds has the advantage that the 
coherence scale is larger than in other mechanisms. As an example, we mention the work 
by Gasperini, Giovannini and Veneziano ^ . which is based on a pre-big-bang cosmology 
inspired by superstrings, which is an alternative to the usual slow roll inflation. The 
dilaton field 4> in the Lagrangian 

J^--V9 ^"^R + \ M d'<P + -/.^F^ (4) 

amplifies the quantum fluctuations of F^^. The magnetic energy spectrum behaves like 
~ kP'^ . The resulting magnetic fields are of the right order of magnitude on a 100 
kpc scale. Recently, however. Turner and Weinberg ^ have argued that this scenario 
requires fine tuning of the initial conditions in order to get enough inflation to solve 
the flatness and horizon problems. 

Bubble Formation at the EW Phase Transition and Magnetic Fields 

In a flrst order EW phase transition bubbles of new vacuum are formed. This 
was used by Baym, Bodecker and McLerran ^ to obtain the generation of a magnetic 
field. The main point is that the bubbles, although overall neutral, have a dipole 
charge layer on the surface, so rotating bubbles generate a field. Although the field 
from each bubble is very small, there is a large number of bubbles, so depending on 
the subsequent development of < S >, in the end a reasonable magnitude can be 
produced. A different mechanism was considered by Kibble and Vilenkin ^: when 
the bubbles collide, the phase of the Higgs field varies, giving rise to currents and a 
magnetic field. Again, in this case one can get a reasonable magnitude provided the 
subsequent development of < S > is favourable. 

Superconducting Cosmic Strings and Other Mechanisms 

There exists a number of other mechanisms for the generation of magnetic fields. 
Vachaspati pointed out that if the gradients of Higgs fields fluctuate, they can induce 
a magnetic field^ at the electrowcak scale. The statistical averaging involved in this 
scenario was discussed in details by Enqvist and me^. As mentioned already in the first 
section, the conclusion is that if line averaging is relevant, one obtains nearly the right 
order of magnitude, since by this mechanism the field at genesis is of order on a 
scale of 1/mw- However, this scenario operates with physically motivated fluctuations 
in gauge dependent quantities like gradients of Higgs fields. Such a procedure is not 
very clear to me. 
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Recently there has been discussions of generation of primordial magnetic fields 
from a network of Witten's superconducting cosmic strings These strings are current 
carrying, and hence produce magnetic fields. It turns out that if the strings are created 
at the GUT phase transition, where the current is very large, they can produce a field 
which is large enough over a sufficient scale, assuming that MHD does not give rise 
to any trouble. On the other hand, superconducting strings created at the EW phase 
transition cannot generate sufficient fields. 

It has also been proposed that a Savvidy-type vacuum, where the energy is lowered 
relative to the trivial vacuum by having a magnetic field , can generate enough field^. 
For SU(A^) the field produced at a temperature T is of order 



At the EW transition (A'^ = 3), this field is far too small. At the GUT transition, 
however, it produces a large enough field for A^ = 5, due to the strong sensitivity of the 
exponent with respect to A^. Whether this field is acceptable depends on the subsequent 
development according to MHD. 

It was proposed long ago by Harrison that magnetic fields could be generated from 
vorticity present in eddies of plasma in the early universe^". This idea was criticized, 
and the eddies were replaced by irrotational density fluctuations by Rees^°. A more 
modern version of this scenario is due to Vachaspati and Vilenkin^°, where the magnetic 
field is generated by vorticity arizing in the wakes of ordinary (i.e. not superconducting) 
cosmic strings. 

Finally we mention that recently Joyce and Shaposhnikov^^ have presented a sce- 
nario which has the potential of leading to quite large fields. The standard model has 
charges with abelian anomaly only (e.g. right-handed electron number) which are essen- 
tially conserved in the very early universe, until a short time before the EW transistion. 
A state with finite chemical potential of such a charge is unstable to the generation of 
hypercharge U(l) fields. Such fields can turn into large magnetic fields, depending on 
their subsequent development. 

It is clear that the physical validity of most, if not all, of these scenarios, depends 
on the subsequent non-linear development of the primordial field, due to MHD. This 
will be discussed in the next section. In the end of this talk we shall also discuss Silk 
diffusion, which is a mechanism for destroying magnetic fields by turning it into heat. 
We shall show that this linear diffusion is in fact counteracted by the non-linear terms 



INVERSE CASCADE FROM MAGNETOHYDRODYNAMICS 

We shall now investigate what happens subsequently to a primordial magnetic field 
generated in the early universe. To simplify things, here we give the details only for 
the non-relativistic case, and mention without giving the arguments, what happens in 
the general relativistic case. 

The Non-Relativistic MHD Equations 

In the rest frame of a plasma consisting of charged particles with current j, we 
have Ohm's law. 




(5) 



of MHD. 



J rest 



a E, 



(6) 
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where a is the conductivity. The universe is a good conductor, so a is very largej^. Thus 
it follows that 

E,est ~ 0. (7) 

In a frame moving with bulk velocity v one therefore has 

E ^ -V X B. (8) 

The induction equation d^/dt = —V x E therefore gives 

dB/dt ^ V X (v X B). (MHD I) (9) 

This is one (out of two) of the fundamental MHD equations. It tells us that the magnetic 
field is influenced by the velocity, and also, if you start from B = no magnetic field 
can be generated. Therefore a seed field is needed in the dynamo mechanism. 

The second fundamental MHD equation is the Navier-Stokes equation with the 
Lorentz force j x B on the right hand side. Here j can be estimated from the Maxwell 
equation j + dE/dt = V x B. The time derivative can be estimated to be small in the 
non-relativistic case|^ so j ~ V x B, so the Navier-Stokes equation with the Lorentz 
force is given by 

dpv/dt + (vV)pv ^ -V(p + ^B^) + (BV)B. (MHD II) (10) 

Here for simplicity we have ignored the viscosity. For a large, this can be generalized 
to the relativistic case at the expense of a considerable increase in the complexity of 
the equations. 

Why and When does MHD have an Inverse Cascade: A Simple Scaling 
Argument 

Now let us suppose that by some particle physics mechanism a primordial magnetic 
field is generated. At the genesis the field has some correlation length, and the crucial 
question is then what happens as time passes. If the correlation length grows smaller, 
corresponding to a cascade, the situation is quite bad, even for the inflationary scenario. 
Such a cascade would appear if the system develops into a more chaotic direction. If, 
on the other hand, we have an inverse cascade, the correlation length grows and the 
system develops towards more order. In an inverse cascade, energy is thus transferred 
from smaller to larger scales. 

It turns out that the situation depends on the initial spectrum. Roughly speak- 
ing, if the spectrum is concentrated at short (large) distances, it will develop into large 
(short) distances. To see this, one can make use of the fact that the MHD equations 
are invariant under the "self-similarity" equations, 

X ^ /x, l^-H, V ^ l\, B ^ /^B, (11) 

■i-In the relativistic era this can be seen from the foUowing estimate: The current is defined by j = rtev. 
The velocity is given essentially by the Newtonian expression v/r w eE/_E, where E is the relativistic 
energy, and r is the average time between collisions. Thus, r w l/na^, where CTx is a typical 
relativistic cross section. Thus, j ~ ne{TeEi/E), so ct w ne^r/E w e^jEoy,. A relativistic cross 
section goes like /T"^ , since the temperature T is a typical momentum transfer. Also, E ^ T . Thus 
(7 « T /e^, which is very large in the early universe, because the temperature is very high. At later 
stages the universe is still a good conductor, for different reasons '^'^ . 

§We have E ^ vB, so dE/dt ^ {v/l)vB ^ {B/l)v'^. But |V x B| B/l, where I is some typical length, 
so the time derivative of the electric field can be ignored relative to the curl of B. 
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and if the viscosity v and Ohmic resistance are included, we further have 

V /i+V, o r^-^o. (12) 

From this it is very easy to show that the magnetic and kinetic energy densities [E-\ 
given in 3+1 dimensions by expressions hke 

£{k, t)V = J d^x £y e'^^^'-y') < B(x, t)B(y, t) >, with V = J £x, (13) 

where 

V J dk £{k,t) ^ y '^^^ < -^^ >^ total magnetic energy, (14) 

must satisfy the scaling relation^^ 

S{k/l,f~H) =1^+^'' S{k,t). (15) 

This is valid in the inertial range, where viscosity and Ohmic resistance can be ignored. 
The general solution of this equation is 

S{k,t) = kP ^{k^^+P^/H), (16) 

with p = —1 — 2h and ^ some function of the single argument k^^'^^^^'^t. The interpre- 
tation of this equation is that, if at the initial time t = the spectrum is k^ (from some 
particle physics mechanism), then at later times it will be governed by the function ip. 
Hence the wave vector scales like 



-2/(3+p) 



(17) 



Thus, if p > —3 there is an inverse cascade, because k moves towards smaller values, 
whereas for p < —3, there is a cascade. Thus, if initially we have a random system 
corresponding to p > 0, then later the system becomes more ordered, as already an- 
nounced. For p = —3 it follows from eq. (15) that the k and t dependence of the energy 
density become uncorrelated. 

In the case when general relativity is included one obtains for a flat, expanding 
universe^^ 

i?(t)^ £{k, t) = F iPik^^+P^^ i), (18) 

where t is the Hubble time, i = J dt/R{t) is the conformal time, t oc \/t, and where 
k is the comoving wave vector, so that the physical wave vector is /cphys = k/R{t). 
Therefore the physical wave vector scales like 

W ~ t-'/('+^V^(i) oc t-(5+f)/2(3+p). (19) 

Thus, if the spectrum starts out with p = 2, corresponding to a Gaussian random 
initial fieldf], we have a scaling of the physical wave vector by instead of t~^'^ 

^If the initial field is given by 

< S,(x,0)Sfe(y,0) >= A (^S.k - ^) S^i^~y), (20) 

then the initial energy is given by 

£{k,0) = X{Su - hki/k^)k^ = 2Afc2. (21) 
Thus the general relativistic scaling g 
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from pure expansion. If p decreases, the effective expansion increases. Tlius, if the 
initial spectrum is characterized by p = 1, we get a scahng fcphys ~ and for p = 

we have fcphys ~ These examples imply physically that the typical size of an eddy 

increases like and t^'^, respectively, on top of the expansion factor. Finally we 

mention that the "scale invariant" initial spectrum with p = —1 {dk/k =scale inv.) 
has an increase of the typical eddy size by t, which means that the eddies follow the 
horizon. 

Numerical Simulations in 2+1 Dimensions 

From the general scaling arguments presented above one cannot deduce the value of 
the scaling function ijj. Here numerical investigations are needed, since realistic analytic 
solutions of MHD are not known. However, a problem arises, since the Reynold number[J] 
is very large in the early universe. For example, in the paper by Brandenburg, Enqvist 
and me^'' the magnetic Reynold number was estimated to be of order 10^^. In numerical 
simulations one cannot reach this value, no matter how much computertime is used. 
We therefore did numerical simulations in 2+1 dimensions with an unrealistically low 
Reynold number^^, taken to be 10. So the non-linear terms are approximately ten 
times as important as the diffusion terms. These terms are, however, needed to achieve 
numerical stability (they act as a short distance cutoff). 

We used the general relativistic MHD equations, which are considerably more 
complicated than the MHD equations discussed in a previous section. We took the 
initial conditions that the velocity vanishes and B is Gaussian random, so the magnetic 
energy spectrum goes like k. Also, we took the energy density to be p = const/i?"^, and 
the pressure p = p/3. 

In fig. 1 the numerical results^^ are displayed. At the initial time there is a rather 
chaotic state, where the magnetic flux lines are either long random walk curves, or 
small closed loops. We used periodic boundary conditions, and satisfied divB = 0. We 
see that in a short time the typical lenght scale increases considerably. In the end of 
the simulation there are quite large eddies. Therefore we clearly see an inverse cascade, 
where order is produced from chaos, in contrast to the usual paradigm. 

Also, the initial velocity v = acquires a spectrum which shows an inverse cascade. 
The velocity is initially induced by the Navier-Stokes equation through the Lorentz 
force. The velocity generated this way then influences the magnetic field through the 
induction equation, etc. etc. 

Numerical Simulations in 3+1 Dimensions: The Shell Model 

As already mentioned, simulations of MHD with large Reynold numbers is not 
possible with present day computers^. The situation gets worse when we go from 2+1 
to 3+1 dimensions. Therefore one needs to make a model which has as many features 
of the real Navier-Stokes as is compatible with practical tractability. In recent years 
the so-called GOY (Gledzer, Ohkitani and Yamada) model has become increasingly 
popular. Another name for this model is the "shell model". It gives results in good 
agreement with experiments, especially as far as the subtle intermittency effects are 
concerned. The model captures a basic feature of turbulence, namely the coupling of 
many different length scales. It is not known whether the model has relation to the real 

"This number can be understood as the ratio between "typical" non-linear terms and the linear vis- 
cosity term. Thus, if Re is large, turbulence is important. 
**This also applies to hydrodynamics, and is perhaps the reason why weather forecasts are pretty bad, 
at least in Denmark. 
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t = A.7 t = 4.7 
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t = 650. t = 653. 
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64x64 128x128 



Figure 1. Left column: magnetic field lines at different times at low resolution (64 x 64 
meshpoints). Right column: magnetic field lines at different times at higher resolution (128 x 128 
meshpoints). This figure is taken from ref. 14. 



Navier-Stokes and MHD. But it nicely illustrates the behaviour of a system in which 
numerical simulations are made difficult by the effect of a huge number of couplings 
between the different length scales. Also, real world conservation laws (energy, helicity) 
are buildt into the model. 

To motivate the model, let us mention that in the Navier-Stokes equations and 
MHD one has terms like 

(vV)v, V X (v X B), (BV)B, (22) 
etc. In Fourier space they e.g. have the form 

(vV)v ^ j d^p Vi{ip){pi - ki)vj{p - k). (23) 

Experience with numerical simulations show that the largest contributions come from 
triangles in k— space with similar side lengths. This is taken as a "phenomenological" 
input in the shell model. 
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At this stage, for numerical purposes, one would discretize k— space. In the shell 
model one of the basic ingredients is a hierachical structure, where |k|— space is divided 
into shells 

kn^y'ko, 71 = 1,2, ...,7V. (24) 

Here A is often taken to be 2. There furthermore exists a complex "velocity mode" 
Vn — v{kn), which can be considered as the Fourier transform of the velocity difference 
|v(|x| + 27r/A) — v(|x|)|. Since /c„ increases exponentially, it covers a wide range of 
corresponding length scales. The model then assumes couplings between neighbours 
and next nearest neighbours, 

(vV)v — > Cij Vn+iknVn+j, (25) 

ij=-2 

where the sum is over neighbours and/or next nearest neighbours to n. The couplings 
Cij in this sum should be made such that energy is conserved in the absence of diffusion. 
Thus, energy conservation 

{v^ + B^)d^x = const (26) 



now corresponds to 
Thus, we need to satisfy 

n=N 



n=N 



J^iKl'^ + W)^ const. (27) 



^ V ^ ^ ^ ^ complex conj . j = 0, (28) 

n=l ^ ^ 

where, as before, t is the conformal time, t = J dt/R(t). In this approach the vectorial 
character is thus lost, but the conservation of energy is kept as an essential feature. 

We should now find equations for the time derivatives respecting the conservation 
of energy. Taking into account some factors from general relativity in an expanding 
universe (the expansion factor as well as the energy density and pressure) we get^^ 

Idvjdi = ikn{A + c)«+,<+2 - 5:+i5:+2) 

+i/c„(i? - iC)«_i<+i - Bl_,Bl^,) (29) 
-ikr,{\B + \A){vl_^vl_, - B*_,B*_,), 

dBjdi = ikn{A - C)«+i5;+2 - 5;+i<+2) 

Mkn{B + \C) (vUB:^, - BUvl^,) (30) 

-ikn{\B - \A){vUBU - ^:-2<-i), 
where with A, B, C arbitrary constants energy is conserved. In 3+1 dimensions, mag- 
netic helicity is also conserved. In the continuum helicity is given by 

H = J d^x AB, (31) 

where A is the vector potential. This conservation is trivial in 2+1 dimensions, since 
there H ~ 0. To mimic conservation of H in the shell model we require that the 
quantity 

n=N 

H,^,,,= Y,{-^TK^BlBr. (32) 

n=l 
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Figure 2. Spectra of the magnetic energy at different times. The straight dotted-dashed Hne gives 
the initial condition (to = 1), the sohd hne gives the final time (t = 3 x fO"'), and the dotted curves 
are for intermediate times (in uniform intervals of A\og{t — io)=0.6). A — 1^ B — —\/2^ and C = 0. 
This figure is from ref. 14. 



is conserved. The reason is that k~^Bn is hke the vector potential. The factor (—1)" is 
a more "phenomenological" factor. The corresponding conservation in hydrodynamics 

(/ ^ v)(i^x =const) has been studied, and it was found that the integrand oscillates 
in sign. This is then taken into account in the shell model by the oscillating factor. 

The requirement that helicity is conserved thus corresponds to taking into ac- 
count 3+1 dimensions, and it leads to the following values for the otherwise arbitrary 
constants A,B,C, 

A=\ B = -1/4, C = 0. (33) 

Using these values, we have 2N coupled set of equations. In our calculations we took 
N —30, corresponding to solving 60 coupled equations. The resulting spectra at differ- 
ent times are shown in fig. 2. Again we see a nice inverse cascade, because as functions 
of the comoving wave vector k the spectra clearly move towards k — 0. 

To give a more precise picture of the change of the spectrum towards large dis- 
tances, we also computed a correlation length defined by averaging over the magnetic 
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10° 10"^ 10^ 10® 10^ 

Figure 3. The correlation lengtli (the "integral scale") as a function of time. The two curves are for 
slightly different models. This figure is from ref. 14. 

energy density, 

(34) 

In turbulence theory this quantity is called the "integral scale" . It is a measure of the 
characteristic size of the largest eddies of turbulence. 

The result is shown in fig. 3. We see that initially the system moves extremely 
rapidly towards larger scales. Clearly MHD (in the shell version) does not like the 
initial Gaussian random state for the magnetic field! The scaling arguments in eq. (18) 
predicts an increase in the eddy size hke This cannot directly be compared to the 
integral scale Iq, since the integrations in eq. (33) are limited by an ultraviolet cutoff, 
which also becomes scaled. However, a fit in ref. 14 gives /q ~ if the steep initial 
increase in fig. 3 is ignored. Taking into account some uncertainty in the fitting, this 
is in good agreement with the scaling in eq. (18). 

EFFECTS OF DIFFUSION: SILK DAMPING 

The effect of diffusion has been ignored in the above discussion, except as a short 
distance cutoff in the numerical calculations. However, this is not realistic, as was 
pointed out by Siegl, Olinto, and Jedamzik^^. This is connected to Silk damping, 
which occurs in the charged plasma because radiation can penetrate the plasma and 
carry away momentum by scattering off the charged particles. Around the time of 
recombination photon diffusion became very important and corresponded to a very 
large photon mean free path. The diffusion coefficient is proportional to the photon 
mean free path, and hence photon diffusion at that time cannot be ignored^^. In a linear 
approximation of MHD it was clearly demonstrated that the magnetic field must be 
destroyed, the magnetic energy beeing turned into heat^^. Silk diffusion would therefore 
remove the hope of understanding primordial magnetic fields from most points of view! 

All hope is not lost, however, since the non-linear inverse cascade, discussed in the 
previous section, counteracts Silk diffusion. While the latter is buzy removing magnetic 
energy at shorter scales, the former is active in removing the energy from these scales 
to large scales. As we have seen in fig. 3, this happens very quickly. Therefore, 
without doing any calculations it is clear that these two mechanisms compete against 
one another. 

To be more precise, numerical simulations arc needed. This was done by Branden- 
burg, Enqvist, and me^^, and the result is that even if Silk diffusion is included, the 



27r 

dk—£(k,t) 
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J dk S{k,t) 
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inverse cascade is strong enough to make the magnetic field survive, at least until close 
to recombination. This should be enough for the dynamo effect to start to operate. We 
refer to the original paper^^ for a full discussion of this. 



CONCLUSIONS 

In conclusion we mention that there are several particle physics models which can 
produce primordial magnetic fields. Of course, they are based on assumptions which 
may not turn out to be ultimately true. For example, there may not be a first order EW 
phase transition, superconducting or ordinary cosmic strings may not exist, etc. etc. 
So when the dust settles, there may not be so many mechanisms which survive. Also, 
it should be remembered that without the inverse cascade, there is no hope to produce 
large enough background fields (this does, of course, not apply to the infiationary 
mechanism), and it may be that for some or all of these models, the inverse cascade is 
not large enough. 

Thus, in estimating the effect of various models one should take into account the 
combined effect of the inverse cascade and Silk diffusion. This will perhaps require 
rather complicated numerical calculations, although some results might conceivably be 
obtained or guessed from simple scaling arguments, as discussed in ref. 13. 

Finally, we mention that there is a very interesting proposal for direct observation 
of a primordial background field^^. The idea is that gamma rays arising from strong 
sources can scatter in a background field, making pair production and delayed pho- 
tons. The spectrum of these photons could then be observed, provided the field is of 
order 10~^^ G or larger^^. If this is technically feasible, important information on the 
spectrum would be obtained, which could then be compared with different models. 
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